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ABSTRACT 

The WMAP satellite, devoted to the observations of the anisotropies of the Cosmic Microwave Back- 
ground (CMB) radiation, has recently provided a determination of the baryonic density of the Universe 
with unprecedented precision. Using this, Big Bang Nucleosynthesis (BBN) calculations predict a pri- 
mordial ''Li abundance which is a factor 2 — 3 higher than that observed in galactic halo dwarf stars. 
It has been argued that this discrepancy could be resolved if the ''Be(d,p)2a reaction rate is around a 
factor of 100 larger than has previously been considered. We have now studied this reaction, for the 
first time at energies appropriate to the Big Bang environment, at the CYCLONE radioactive beam 
facility at Louvain-la-Neuve. The cross section was found to be a factor of 10 smaller than derived 
from earlier measurements. It is concluded therefore that nuclear uncertainties cannot explain the dis- 
crepancy between observed and predicted primordial ^Li abundances, and an alternative astrophysical 
solution must be investigated. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — stars : Population II — cosmo- 
logical parameters, early universe 



1. INTRODUCTION 

Us ing the WMAP -determination of th e bary onic den- 
sity lIBennett et a,1. \VEM ISnergel et al. II200I. one ob- 
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tains predictions of the abundances of the light ele- 
ment isotopes pr oduce d in Big Bang Nuc leosynthesis 
jCoc et al. 1127^1 l20nl ICvburt et al. II2003D . While the 
overall values from theoretical predictions and from the 
observational determinations of the abundances of D and 
*He are in good agreement, the theory tends to predict a 
higher ^Li abundance (by a factor 2 to 3) tha n is observed 
in the atmospheres of halo dwarf stars l|R,van et al. I 
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l2000(l . The NACRE compilation ijAngulo et al. llT999(l 
provided a new set of reaction rates that were used to 
update the predictions of contemporary Big Bang Nucle- 
osynthesis (BBN) ( Vaneioni-Flam ct al. 2000). At that 
time, the baryonic densities obtained from CMB observa- 
tions on the one hand and comparison between BBN cal- 
culations and spectroscopic data on the other hand were 
only marginally compatible llCoc et a l. ||2002 |). In order 
to im prove the nuclear network. IDescouvemont et al. I 
l)2004tl recently performed a re-analysis of low energy 
data from the 10 key nuclear reactions involved in 
BBN, by using R-matrix theory l)Lane fc Thomas1ll95S() 
and evaluating the remaining uncertainties in a statis- 
tically robust for malism. Using this improved network, 
iCoc et all l)2004j) have recently calculated BBN light ele- 
ment productions assuming for the baryonic density the 
very precise value provided by WMAP ( So ereel et al. I 
l2003|l and obtained ^Li/H= 4.15l:|^;4^ x 10-^° compared 
to the observed value Li/H ~ 1 — 2 x 10^^", confirming 
the ^Li discrepancy. 

However, it has been shown l|Coc et al. II2004I) that the 
''Be(d,p)2a reaction (which destroys the ^Be that is the 
source of ^Li at high baryonic density), would solve the 
''Li problem i/its cross section were much higher than as- 
sumed. Importantly, prior to the present work, no direct 
experimental data at BBN energies were available (for 
T = 0.5 - 1 GK, the Gamow window is = 0.11 - 0.56 
MeV). In fact, the ^Be(d,p)2a reaction rate relied on an 
extrapolation made by Parker ( 1972) based on experi- 
mental data at center-of-mass (c m.) energies of 0.6 to 
1.3 MeV from iKavanagh I l)196CI|) . In this experiment, 
protons corresponding to the ^Be 0+ ground state (g.s.) 
and first excited state (3.03 MeV, 2+) were detected at 
90° using a Nal(Tl) detecto r . Assu ming an isotropic an- 
gular distribution, iParker I l)1972f) multiplied the mea- 
sured differential cross section by An and by a further 
factor of 3 to take into account the estimated contribu- 
tion of the higher energy ^Be states, not observed by 
^avanagh (1960). Consequently, a constant 5- factor of 
100 MeV-barn was adopted. 

In order to obtain ^Be(d,p)2a reaction cross section 
at BBN energies, we have performed an experiment at 
the CYCLONE radioactive beam facility at Louvain-la- 
Neuve, Belgium, using an isobarically pure ''Be radioac- 
tive beam. The experimental method and results are 
presented in Section 2. The astrophysical consequences 
are discussed in Section 3. The conclusions are given in 
Section 4. 

2. EXPERIMENTAL METHOD AND RESULTS 



Fig. 1. — Schematic view of the experimental set-up. 



The measurements were performed using a post- 
accelerated ''Be^+ radioactive beam at a nominal en- 
ergy of 5.8 MeV provided by the CYCLONEllO cy- 
clotron. A detailed description of the production of the 
''Be beam can be found inlGaclcns ct al. ( 2003). To sup- 
press the contamination from the '^Li isobaric beam, the 
''Be beam was completely stripped to ''Be'*+ by trans- 
mission through a thin ^^C foil, prior to analysis by a 



dipole magnet. Prior to the ''Be(d,p)2a measurement, 
the beam energy was determined using a calibrated Si 
detector situated at 0°. A laboratory energy of 5.55 MeV 
(FWHM ^ 4%) was determined, including a correction 
for pulse height defect. This energy was degraded to 
1.71 MeV (FWHM ~ 12%) using a 6 fim. Mylar foil lo- 
cated at 50 cm upstream of the target. No ''Li contam- 
ination was observed, consistent with lithium isotopes 
being unable to support a 4+ charge state. The target 
consisted of a 200 /ig/cm^ (CD2)„ self-supporting foil. 
With this set-up, we were able to investigate the center- 
of-mass energy range between 1.00 and 1.23 MeV (for 
a beam energy of 5.55 MeV, without degrader) and be- 
tween 0.13 and 0.38 MeV (for 1.71 MeV, with degrader). 
The cross section measurement was averaged over these 
energy ranges. In addition to the feeding of the ground 
and first excited states of *Be l)Kavanagh I l960'l. we were 
able to observe the ''Be-l-d reaction via other kinemati- 
cally allowed higher energy levels, mainly through a very 
broad 4+ state (F ~ 3.5 MeV) situated at an excitation 
energy of 11.4 MeV in ^Be (Tillev et al. 2004). At the 
beam energy of 5.55 MeV, several states in *Be above the 
''Be-l-d-p threshold are present but due to the Coulomb 
barrier in the final state, their contribution are expected 
to be negligible. The Q value of the ''Be(d,p)^Be reaction 
is 16.49 MeV, thus the laboratory energies of protons and 
a particles are high. For example, a 5.55 MeV ''Be beam 
traversing a 200 ^g/cm^ (CD2)„ target will lead to the 
production of protons with energies anywhere between 
about 7.5 and 22 MeV, for the range of angles covered. 
Thus, to clearly identify the protons from the ''Be-l-d re- 
action from those arising from reactions on the C content 
of the targ et, a stack of two 'LEDA' silicon strip detec- 
tor arrays l|Davinson et al. I l2000) were employed cover- 
ing a laboratory angular range of 6'iab = 7.6° — 17.4°. 
A schematic view of the experimental set-up is shown 
in Figure 1. The AEi detector consisted of eight sec- 
tors of 0.3 mm thickness, while the AE2 detector con- 
sisted of four sectors of 0.3 mm thickness and four of 
0.5 mm thickness. They were calibrated using a 3-line 
a-source (^^^Pu, ^^^Am, ^^''Cm) and a precision pulser. 
This AEi — AE2 detector system allowed a clear identifi- 
cation of the protons produced in the ''Be(d,p)^Be reac- 
tion. All the particles that are not stopped in the front 
AEi detector and that are either stopped or that left 
energy on the back AE2 detector are protons (having an 
energy of more than 6.5 MeV) that populated levels up 
to the 11.4 MeV state in ^Be. We were able to measure 
the ''Be(d,p)2a cross section up to an excitation energy 
in ^Be of E^ = 13.8 MeV for a beam energy of 5.55 MeV 
and of E^ = 11.5 MeV for 1.71 MeV. Only about 50% 
of the contribution of the 4+ broad state was observed 
at 1.71 MeV. Other light particles (p, d, ^He, '^He) from 
''Be-|--'^^C reactions, as well as recoils and scattered par- 
ticles, were completely stopped in AE'i. 

Figures 2 and 3 show spectra obtained at beam ener- 
gies of 5.55 and 1.71 MeV, respectively. The spectrum 
obtained at 1.71 MeV (Figure 3) was accumulated over 
about 26 hours of running time with an averag ed ''Be 
beam intensity of 2 x 10^ pps. As can be seen, the 
proton signals are well separated from the background 
signals {AE2 < 1 MeV), which are produced by ran- 
dom coincidences of a particles, scattered ''Be and re- 
coil ions in the AEi detector with (3 particles in Ai?2. 



The locus with negative slope contains protons that have 
passed through the front AEi detector and stopped in 
the back AE2 detector. The two loci with positive slope 
are events in which the proton has sufficient energy to 
pass through both detectors completely. There are two 
bands because of the two different thickness AE2 detec- 
tors. The most strongly populated regions at the lower 
left of these bands correspond to protons losing the least 
energy in passing through the silicon, and thus to the 
highest energy events. By considering the kinematics 
and energy losses in silicon (^SRIM 20 03), together with 
the straggling of the beam and experimental energy reso- 
lution, one may then identify events on the positive slope 
locus up to 2.5 (3.9) MeV for the 0.3 (0.5) mm Si thick- 
ness wafer, as corresponding to events in which the re- 
coiling ^Be nucleus is in either the 0+ g.s. or the 2+ 
excited state (the energy resolution is insufficient to re- 
solve the two). The total statistical error was 10% for the 
beam energy of 1.71 MeV and less than 2% for 5.55 MeV 
(for protons populating the 0+ and 2+ states it was 13% 
and 2.5%, respectively). The absolute normalization was 
obtained using events arising from the elastic scattering 
of the ^Be on the C content of the target (as recorded 
by the AEi detector in which the ''Be are stopped), and 
assuming that the ^Be-|-^^C elastic scattering follows the 
Rutherford law. This assumption is realistic at energies 
below the Coulomb barrier, as is the case here. 

To calculate the average cross section over the energy 
ranges and angular coverages (da/dft), the number of 
counts was corrected for the detector solid angle (un- 
certainty ±5%), the number of deuterons in the target 
(±10%) and the total number of incoming beam particles 
(±7% and ±26% at the higher and lower beam energies, 
respectively), and transformed into the cm. system. 




Fig. 2. — AE1-AE2 spectrum at a beam energy of 5.55 MeV on 
a 200 ^gr/cm^ (CD2)n target. The cm. energy range covered is 
1.0 to 1.23 MeV. 

The proton angular distribution over the angular range 
covered here was found to be isotropic at both energies. 
Thus, we assumed full isotropy and calculated the aver- 
age cross sections, ct = 7.5 ± 0.8{stat) ± 2.6{sys) mb at 
the effective energy of 0.37 MeV and ct = 386 ± 7(stat) ± 
50{sys) mb at 1.15 MeV. The summed contribution of 
the 0+ and 2+ states was about 64% of ct at 1.15 MeV. 
At 0.37 MeV, a includes the contribution of the ground 
and 2+ states and about 50% of that of the 4+ broad 
state. Due to the low penetration probability (£ = 4), 
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Fig. 3. — Same as Fig. 3 for a beam energy of 1.71 MeV, corre- 
sponding to a cm. energy range of 0.13 to 0.38 MeV. 



the contribution of the 4+ state should be less than 36%. 
Thus, the a value at 0.37 MeV corresponded to more 
than 80% of the total cross section (for a 4"*" state with 
r ~ 3.5 MeV). This was taken into account in the sys- 
tematic uncertainty. 

In nuclear astrophysics it is usual to present the cross 
section in the form of the astrophysical S'-factor S{E) 
given by IfClavt on 1983), 

S{E) ^ a{E) exp(2nri)E, (1) 

where 77 is the Sommerfeld parameter (77 = ZiZ2e^ /hv, 
with Zi and Z2 the charge numbers of the target and 
beam and v is the velocity) and E is the effective 
cm. energy. In the absence of sharp resonances, the S'- 
factor varies smoothly with energy. Figure 4 shows the 
''Be(d,p)®Be astrophysical S'-factor S{E) in MeV-barn as 
a funct ion of the cm. ene rgy. For a comparison with the 
data of lKavanagh I (|196(l ) (open circles) , the present data 
including only contributions from the ground and first ex- 
cited states of ®Be (filled circles) are shown. The agree- 
ment with thepK^avanash ( 1960) data at overlapping en- 
ergies is satisfactory, given the systematic uncertainties. 
The total S-factor is also shown (filled triangles). The 
present d ata show that the higher energy states not ob- 
served bv lKavanagh 1 1)1960|) contribute about 35 % of the 
total S'-factor instead of the 300% estimated by IParker I 
ljl972() . Hence, the ^Be(d,p)®Be reaction rate is smaller 
by a factor of about 2 at energies in the range 1.0 to 
1.23 MeV and by about 10 at energies relevant to BBN, 
than previously estimated. This excludes a nuclear so- 
lution to the primordial lithium abundance problem via 
the ^Be(d,p)*Be reaction as its effect is completely neg- 
ligible compared to the 7% (Icr) nuclear uncertainty on 
the ''Li yield. Nevertheless, these results allow a more 
accurate determination of the ^Li abundance using BBN 
models. 

3. ASTROPHYSICAL CONSEQUENCES 

Since the pioneering work of ISnite S^^ ((1982), 

who found a value of Li/H« 1.2 x 10~^° independent 
of Fe/H (for [Fe/II]< —1.3) there have been many inde- 
pendent observations of Li confirming the existence of 
a plateau and suggesting that this abundance reflects 
the primordial Li value. However, the Li abundance 
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Fig. 4. — Astrophysic al S factor of the ''Be(d,p)*Be reaction. 
Open circles: data from i Kavanaghl 119601 : filled circles: present 
data including contributions from the ground and first excited 
states of *Be only; filled triangles: total 5-factor derived from the 
present experiment. The vertical error bars are the total error. The 
horizontal dotted bars indicate the energy range covered at each 
data point. The Gamow window for a typical BEN temperature 
T=0.8 GK is also shown. 



extracted from observations depends dras tically on the 
assum ed surface temperature of the star (iFields et al. I 
1)200,^ V Recent observations l)R van et al. I l200fl() have 
lead to Li/H (1.23;^-^ ^) x IQ-^" which is very 

[^82). The 



Li/H (1.23^ 
close to the first evaluation 



iiie IV, r)Dr 

more recent work studied and quantified the various 
sources of uncertainty: extrapolation, stellar depletion 
and stellar atmosphere parameters. Compared to the 
WMAP+BBN value, the discrepancy is a factor of ^ 3.4. 
If it is shown that there is a mechanism by which the 
outer layers of Population II stars are transported deep 
into the stellar interior, then there are several ways in 
which Li abundances might be depleted over the life- 
time of the star. In this context, the current esti- 
mates for possible depl etion factors m ay be i n the range 
~ 0.2-0.4 dex ( Pinson neault et al. 1 2002: Richard et al. I 
120041 iVauclaiT fc Char bonnel 1 11998). However, the data 
typically show negligible intrinsic spread in the Li abun- 
dance leading to the conclusion that depletion in these 
stars is of the ord er of 0.1 dex. 

Recently, iMelendez fc Ramirez I ()2004() have obtained 
a higher value for the Li plateau abundance (2.34x 10^^°) 



due to a new effective temperature scale which is higher 
at low metallicity. This new evaluation diminishes the 
discrepancy, without canceling it. The observation of 
^Li is also of interest, since, because it is more frag- 
ile than ^Li, it can provide yet more s evere constraints 
upon possible deple tion mechanisms ijLambert 1 120041: 
iRollinde et 'ani2005j) . Finally, in spite of the various 
uncertainties related to Li observations and to the stel- 
lar models, it is very difficult to reconcile the BBN ^Li 
and the Spite plateau which presents a narrow dispersion 
all along the metallicity scale. 



4. CONCLUSIONS 

The existence of the Spite plateau for Li seems to indi- 
cate that low metallicity halo stars are indeed representa- 
tive of the primordial BBN abundance. In particular, the 
isotope ^Li plays a key role as a bridge between Big Bang 
Nucleosynthesis, stellar evolution and galactic cosmic- 
ray nucleosynthesis. At present there is a significant 
discrepancy between the BBN-predicted ^Li abundance 
(assuming a baryon density consistent with the concor- 
dance model derived from observations of anisotropics 
in the microwave background) and the abundance deter- 
mined from the observations of Li in the atmospheres of 
halo stars. The experiment reported here demonstrates 
that the ^Be(d,p)2 a iS-factor at BBN energies was not 
underestimated bv iParkerl l)1972|) but, on the contrary, 
overestimated. The discrepancy cannot therefore be re- 
solved by nuclear physics inputs to BBN calculations. 
The remaining conventional options (those not invoking 
physics beyond the Standard Model) are an adjustment 
of the stellar input parameters needed to extract the Li 
abundances from observations, or stellar depletion of ^Li. 
However, models must be constructed to avoid disper- 
sion in the ^Li abundances over a wide range of stellar 
parameters, which is a real challenge. The origin of the 
discrepancy in the Li abundance remains a challenging 
issue. 
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